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Diffusivity is a key quantity in describing velocity fluctuations in granular materials. These 
fluctuations are the basis of many thermodynamic and hydrodynamic models which aim to provide 
a statistical description of granular systems. We present experimental results on diffusivity in dense, 
granular shear flows in a 2D Couette geometry. We find that self-diffusivities D are proportional to 
the local shear rate 7 with diffusivities along the direction of the mean flow approximately twice as 
large as those in the perpendicular direction. The magnitude of the diffusivity is D ~ 7a 2 where a is 
the particle radius. However, the gradient in shear rate, coupling to the mean flow, and strong drag at 
the moving boundary lead to particle displacements that can appear subdiffusive or superdiffusive. 
In particular, diffusion appears to be superdiffusive along the mean flow direction due to Taylor 
dispersion effects and subdiffusive along the perpendicular direction due to the gradient in shear rate. 
The anisotropic force network leads to an additional anisotropy in the diffusivity that is a property 
of dense systems and has no obvious analog in rapid flows. Specifically, the diffusivity is supressed 
along the direction of the strong force network. A simple random walk simulation reproduces the 
key features of the data, such as the apparent superdiffusive and subdiffusive behavior arising from 
the mean velocity field, confirming the underlying diffusive motion. The additional anisotropy is not 
observed in the simulation since the strong force network is not included. Examples of correlated 
motion, such as transient vortices, and Levy flights are also observed. Although correlated motion 
creates velocity fields which are qualitatively different from collisional Brownian motion and can 
introduce non-diffusive effects, on average the system appears simply diffusive. 



I. INTRODUCTION 



Overview 



Despite the prevalence of granular materials in nature 
and industry, a coherent understanding of granular flows 
is still lacking. Particularly in dense systems, features 
such as jamming, shear bands, and the coexistence of 
solid- and liquid-like regions make it difficult to offer 
a simple theoretical description. Fluctuations in both 
the force network and particle velocities can be of the 
same magnitude as the mean values and are known to be 
important aspects of the microscopic behavior of dense 
granular flowspj . Due to the complexity of these systems, 
one of the key goals of current research is to develop a 
statistical description of steady state behavior, such as a 
thermodynamic or hydrodynamic model. Fundamental 
to statistical approaches is understanding the mean fluc- 
tuating part of particle motion, which is described by a 
granular diffusivity. 

Sheared granular systems have received considerable 
attention recently [lj, |2j, |3j as an important example of 
granular flow. Diffusion, in particular, has been studied 
in a variety of granular systems, such as vibrated grains 
HIH, tumblers HH, chute flows [ 111 [Toll , and sheared 
systems Il8lll9f . but these stud- 

ies have predominantly focused on rapid flow regimes. 
Understanding slow, high density, flow is not trivial 0] ; 
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there is no replacement at the fundamental level for colli- 
sionally based kinetic theories that are expected to apply 
only in the dilute rapid flow regime. 

In this paper, we characterize the self-diffusivity of 
grains in a 2D Couette shearing experiment by study- 
ing individual particle trajectories over time. In contrast 
to most previous results, we focus on quasistatic dense 
flows. In this regime, since particles are constantly in 
contact with their neighbors, interactions are not colli- 
sional and material flow is largely confined to a shear 
band with a nominal thickness on the order of 5 particle 
diameters. 

Several observations from the present experiments are 
noteworthy: 1) We find that particle diffusivity is pro- 
portional to the local shear rate, with diffusivities ap- 
proximately twice as large along the mean flow direction 
as the perpendicular direction. 2) We show that unlike 
rapid flows, the anisotropic force network induces a sub- 
stantial anisotropy in the diffusivity. This is in addition 
to the usual anisotropy induced by the direction of mean 
flow. 3) Care must be taken when calculating diffusiv- 
ities in a shear gradient |20| . as motion can appear to 
be subdiffusive or superdiffusive due to a gradient in the 
shear rate or Taylor dispersion |2l| . We show through a 
simple Fokker-Planck model that apparent sub- or super- 
diffusive behavior can be attributed to shear gradient and 
boundary effects. 4) Hence, the grain motion is statisti- 
cally consistent with a simple random walk in the pres- 
ence of shear gradients. 5) Nevertheless, at larger spa- 
tial scales, we occasionally observe correlated motion and 
Levy flights. But these events are rare and do not have 
a significant impact in the mean. 
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B. Models for Granular Diffusion 

In the kinetic theory approach, a granular temperature 
is often defined as T oc ({v— v) 2 ), with instantaneous par- 
ticle velocity v and mean velocity v, in which the velocity 
fluctuations contribute to a temperature in analogy with 
molecular gases. 

A different approach was recently proposed by Makse 
and Kurchan, who applied uniform shear in a numerical 
experiment and measured diffusivity D and mobility \ 
to define a temperature D/x by analogy with fluid sys- 
tems [l4|. I n their model, they report that the th law 
(thermal equilibration) is satisfied in a bidisperse mix- 
ture, supporting the thermodynamic picture. This is in 
contrast to experimental measurements of kinetic granu- 
lar temperature which find a lack of e quip artition when 
different types of particles are present [22, 23]. 

Isotropic Brownian diffusion in an unbounded system 
is often characterized by the time evolution of the second 
moments of a probability distribution function (PDF). 
For example, 



as above, and a diffusion tensor, D, with elements D x 
D yy and D xy — D yx where 



(x 2 ) = 2Dt 



(1) 



where x is the particle position relative to its initial po- 
sition (v = Ax /At for a small time step At), D is the 
diffusivity, and t is time. 

More generally, diffusion must be described by a ten- 
sor. For instance, diffusivities along the flow direction 
in granular gase s are in general different from transverse 
diffusivities (121 Il7l Il8| . 

Diffusion in even a simple shear flow is complicated 
by Taylor dispersion effects pf|. in which diffusive mo- 
tion couples to the mean flow leading to larger dispersion 
along the flow direction, as recently elucidated in systems 
of noncolloidal particles 20]. In this case, (x 2 ) is nonlin- 
ear in time, i.e. it contains higher order corrections due 
to the coupling of the shear to the diffusive motion. Sim- 
ply subtracting the mean flow from particle trajectories 
and computing diffusivities does not give accurate results 
in this system |20|. 

In particular, for flow of the form v — jyx, i.e. uniform 
unbounded shear flow in two dimensions in which there 
is a constant shear rate 7 creating a velocity gradient in 
the y d irection, the second order moments are given by 

Mm 



2Dyyt 



(yy) 

(Xy) = 2D X yt + Dyyjt 2 



(2) 
(3) 



(xx) 



2D XX t + 2D X yjt 2 + -Dyyj H (4) 



These equations describe an ensemble average of particle 
positions (x(t),y(t)) relative to the particle's initial loca- 
tion (i.e. ((x(0), y(0)) = (0,0)), without subtracting the 
mean flow. 

These relations follow naturally for a PDF described 
by Brownian anisotropic diffusion with mean local flow v 



dP/dt — v ■ VP + V • DVP. 



(5) 



Here, x corresponds to the streamwise direction and 
y is perpendicular to x. The left hand sides of Eq. 12141 
correspond to the y 2 , xy, and x 2 moments of P. 

The higher order terms in Eq. ^ are due to Taylor dis- 
persion. For instance, the i 3 term arises because diffu- 
sive motion along ±y moves grains to regions of different 
mean velocity v(y)x which tends to increase their sepa- 
ration or dispersion along the x direction. These higher 
order terms contribute to mean squared displacements 
which therefore appear superdiffusive. 

It must be emphasized that equations 12141 are derived 
assuming that the diffusivities and shear rate are con- 
stant in space and time over an infinite plane. This con- 
dition is not met in the current experiment which changes 
the specific form of the correction terms. 

Several issues concerning a diffusive picture must be 
addressed for sheared dense granular materials to deter- 
mine whether Brownian diffusion applies. Two of these 
issues are the presence of a shear band and the limiting 
boundary at the shearing surface. Even assuming that a 
diffusive description is applicable, it remains to be deter- 
mined on what temporal or spatial scales such a descrip- 
tion should apply. In dense quasistatic flows, grains are 
generally close to a jammed state in which particles are in 
constant contact. Motion of grains requires the creation 
of voids, so correlated motion might be expected to be 
particularly important in dense 2D systems where paths 
are constrained. Short-lived vortex structures have been 
seen in 2D granular simulations and experiments on 
2D shearing of foams[25|. These are potential deviations 
fom Brownian diffusive behavior which might affect the 
time evolution of the moments. 



C. Previous Measurements and Simulations 

The full diffusion tensor has been measured in granular 
gases usin g kin etic theory |l7j , simulations of rapid gran- 
ular shear , and shearing of non-colloidal suspensions 

Substantial work on granular diffusivity in rapid flows 
has been done by Hsiau and coworkers who have mea- 
sured self-diffusion coefficients in a variety of granular 
systems [U EE3, Ez| . They find that fluctuations are 
anisotropic, with the largest fluctuations along the flow. 
Diffusivities are found to increase with shear rate and 
depend on the square root of the granular temperature 
T in agreement with kinetic gas theory. Other results 
in a similar chute flow were subsequently presented by 
Natarayan et al.\ld^. 

Losert et al. studied a 3D fluidized Couette experi- 
ment , in which velocity fluctuations were found to be 
slightly larger in the direction along the mean flow. These 
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fluctuations decrease roughly exponentially far from the 
inner cylinder, but decrease more slowly than the average 
velocity. 

Radjai and Roux studied particle velocity fluctuations 
in numerical simulations under homogeneous strain in 
which there was no shear band formation [15| • They mea- 
sured anomalous diffusion with an exponent of 0.9 (rather 
than 0.5 for ordinary diffusion) which they attributed to 
long-time configurational memory of a granular medium 
in quasistatic flows. 

Diffusivities have also been measured in a 3D rotating 
tumbler 0,0, 2D swirling flow 0, chute flow 0, sim- 
ulations of shaken spheres p|, and simulations of small 
numbers of spheres in suspension |29| . Earlier studies 
primarily addressed rapid flows from kinetic theory |16| . 

Although these studies are relevant here, we note that 
the displacements were assumed to be purely induced by 
the shear flow and no attempt was made to investigate 
the role of the force chain network. 



D. Organization of Presentation 

The paper is organized as follows. In Section^ we de- 
scribe the experimental techniques. We present diffusion 
measurements in Section IIIII and results from a random 
walk simulation in Section llVl We show the impact of 
the anisotropic force network in Section We discuss 
diffusivities determined from the velocity autocorrelation 
functions in Section IVll In Section fVIII we show exam- 
ples of intermittent vortices and Levy flight trajectories, 
and in Section IVIIII we draw conclusions. 



II. EXPERIMENTAL TECHNIQUES 

The experiment is performed with a 2D Couette appa- 
ratus, as sketched in a top view in Fig. ^ The granular 
material (B) consists of a bidisperse mixture of about 
40,000 disks (diameters d s = 0.42 cm, d L = 0.50 cm, 
thickness = 0.32 cm) in a ratio of 3 small : 1 large. The 
bidisperse mixture is used to inhibit crystalline order- 
ing of the disks. The disks lie flat on a Plexiglas sheet 
bounded by an outer ring (i? Q = 51 cm)(C) and an in- 
ner shearing wheel (Ri = 20.5 cm) (A). A Plexiglas sheet 
covers the experiment to protect the experiment from ex- 
ternal perturbations, but the sheet does not contact the 
particles. The shearing wheel is rotated at a frequency 
/ of 0.1-10.0 mHz or a speed of v « 0.013 - 1.3 cm/s 
at the shearing surface. The experiment is initially run 
for at least one revolution of the shearing wheel in order 
to avoid effects from transients, an issue that will be ad- 
dressed in another paper |30| . The shearing wheel and 
the outer ring have teeth with gaps comparable to the 
size of the smaller particles. 

The system is lit from below and observed from above 
using a 2 megapixel CCD camera at a frame rate of up 
to 7.5 Hz. Sequences of approximately 1500 images are 




3.9 cm 



102 cm 



FIG. 1: Schematic of experiment as viewed from above. The 
granular material (B) is contained by the shearing wheel (A) 
and the outer ring (C). Their radii are Ri = 20.5 cm and 
Ro = 51 cm respectively. On the right is a section of an 
experimental image of the grains. 

analyzed to determine particle trajectories. Disks within 
about 20 particles diameters of the shearing wheel are 
marked with lines, thus allowing us to track particle po- 
sition and orientation and to identify particles by size. 
Approximately 4000 grains are typically in the field of 
view, of which 2500 are marked. Images are threshholded 
and the orientation and position of lines on the disks are 
found. Sufficient temporal resolution is used such that 
each particle in a frame can be connected with the clos- 
est grain in the subsequent frame to establish particle 
trajectories. 

The force network can also be visualized since the 
grains are made of a photoelastic material |3l|. When 
polarized light travels through the disks, it experiences a 
phase shift (birefringence) proportional to the difference 
in principle stresses <t 2 — When the disks are illu- 
minated between crossed polarizers, grains under larger 
stress are seen as regions of larger gradients in light in- 
tensity. In this way, the force network is visualized as a 
network of bright lines on a dark background. Additional 
details were presented by Howell et al.[j^. The polariz- 
ers are removed for measuring particle trajectories and 
diffusivities. 



III. DIFFUSION MEASUREMENTS 

A. Mean Velocity Profiles 

We first consider the mean properties of the flow. In 
Fig. |21 we show the mean tangential velocity, v§ , versus 
radial distance from the shearing surface, r = R — Ri, 
where R is the distance from the center of the shear- 
ing wheel. The velocities are scaled by the velocity of 
the shearing surface Vb = 0.28d/s where d is the mean 
particle diameter (d — (ds + di,)/2). For this particu- 
lar run, we used a frame rate of 1.08 Hz, so that the 
shearing wheel was displaced 0.25 d between each of the 
1080 images, which, in total, correspond to one revolu- 
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V = 0.925 exp(-0.284 r - 0.0534 r 

- - V = 1.071 exp(-0.521 r) 
a experiment 
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FIG. 2: Mean tangential velocity vs. radial distance from the 
shearing surface plotted on a (A) linear and (B) logarithmic 
scale for a particular run. The velocity v is scaled by the 
velocity of the shearing surface Vo- The experimental data 
(A, r < 7.5 d) is fitted to an exponential (dashed line) and a 
gaussian (solid line). / = 1 mHz and Vo = 0.28d/s. 
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FIG. 3: (A) Tangential velocity versus radial distance from 
the shearing surface using experimental runs at different 
frame rates (= i) as described in text. Velocities v are 
shown relative to the the velocity of the shearing surface 
Vo = 0.28 d/s. (B) Data is shown where velocities can be 
resolved given the frame rate and number of pictures. An 
exponential tail is observed in which ve oc e _1 ' 25r . 



tion of the shearing wheel. The limiting value for vg of 
approximately 10~ 4 corresponds to the sensitivity of the 
measurement for a typical number of images and spatial 
resolution, e.g. for an image resolution of 20 pixels per 
diameter, a grain displacement of 1 pixel over the entire 
run would give a mean velocity vg = ■ ^ • 1.08 Hz 
= 5 x 10~ 5 ^ (or v/Vq m 10~ 4 ) and velocities smaller 
than this cannot be resolved. Motivated by previous re- 
sults for Couette shear we fit the data (r < 7.5 d) 
to an exponential (vg(r) = 1.071 exp(-0.521 r)) and to a 
gaussian (v e (r) = 0.925 exp(-0.284 r - 0.0534 r 2 )) . 

From this individual run, the Gaussian fit seems most 
appropriate. However, this is an artifact of the data reso- 
lution. Additional data with slower frame rates points to 
an important issue concerning particle tracking velocime- 
try. If the number of images N remains fixed, by taking 
data using slower frame rates, velocities of slower parti- 
cles further from the shearing surface can be accurately 
measured while faster particles at the shearing surface 



can no longer be accurately tracked. We show these re- 
sults in Fig. |3 in which the frame rate (= 1/At) is varied 
for different runs at the same imposed shear rate (/ = 1 
mHz, N — 1500). The individual curves are accurate over 
a particular range of velocities based on the frame rate 
and number of pictures in the run. In Fig.[3}3, we show 
data for each set within this range. It becomes evident 
that the velocity profile has an exponential tail which is 
obscured when simply analyzing a single run. Previous 
results have shown exponential [lj, [l9j , Gaussian |2| , and 
similar strongly decaying || velocity profiles for Cou- 
ette flow. Authors of these studies 0, |2] have suggested 
that the differences in measured profiles may depend on 
whether the flow is 2 or 3D or on whether the parti- 
cles are rough or smooth. The present data suggest that 
an additional factor may be spatio-temporal resolution. 
Particle tracking issues in particular have been addressed 
recently by Xu et aL[33j. 

We conclude that correct tracking occurs for velocities 



that approximately satisfy 
1 1 d 



(6) 



where N = number of pictures, and the image resolution 
is 20 pixels per diameter d. The lower speed limit is set 
by image and temporal resolution. The upper limit is 
chosen to resolve the occasionally fast displacements well 
above the mean. Note, however, that the upper cutoff 
is not an issue when the speed of the shearing surface is 
less than d/ At, since generally, all particle displacements 
can be resolved, e.g. for At = 0.3 s and 0.925 s, the 
requirement that Ug < 0.1 is not necessary. 
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B. Radial and Tangential Diffusivities, D rr and Dgg 

We measure diffusivities by tracking individual parti- 
cles, and hence their displacements, over time in both the 
radial (f) and tangential {&) directions. We subtract the 
mean flow (Fig. [5J from the tangential velocity compo- 
nent at each time step. 

Although Eqs. 1 2141 characterize absolute x and y dis- 
placements without subtracting the mean flow, they are 
also predicated on a velocity profile such that (x(t)) = 
(y{t)) = which is true for a uniform shear rate in an 
infinite domain, neither of which are true here. The ex- 
ponential velocity profile observed in the experiment in- 
validates this and prevents us from directly comparing to 
Ecis. 12141 Deriving the corresponding moment evolution 
equations with the exponential profile of the velocity and 
diffusion fields is significantly more complicated and be- 
yond the scope of this paper. Instead, we use a random 
walk model in Section llVl to model the system. Nonethe- 
less, the theory above explains the origin of the Taylor 
dispersion effects that we observe. In order to avoid the 
effects of a mean displacement due to the locally varying 
shear rate, we remove the mean flow. The resulting mean 
displacement squared is plotted versus time and averaged 
for different particles initially within the same radial bin 
(bin size = d or d/2). An initially linear evolution in- 
dicates ordinary diffusive behavior with the slope of the 
line equal to 2D. 

Fig. 0] shows a typical example of the mean displace- 
ment squared for the tangential ((RA9) 2 ) and radial 
((Ar) 2 } directions for particles in the shear band. Here, 
the A notation reminds us that the mean flow is sub- 
tracted from the data. The dotted lines are linear fits 
for t < 30s giving diffusivities proportional to the slopes. 
The tangential diffusivity is approximately double the ra- 
dial diffusivity at small times. The former is expected to 
deviate from a straight line due to the higher order terms 
similar to those in Eq.0J Note, however, that for small t, 
the linear term in Eq. ^dominates, and Taylor dispersion 
effects are not present. 

One might worry that using early times would be inac- 
curate when diffusivity is generally defined as a long time 
behavior. In particular, results for more rapid flows show 



FIG. 4: Mean displacement squared vs. time for tangential 
and radial directions for particle trajectories starting at 2d < 
r < 3d and shearing wheel frequency of / = lmHz. Dotted 
fir 




Time (s) 

FIG. 5: Velocity autocorrelation versus time for particles in 
2d < r < 3d for radial and tangential velocity components. 



an initial ballistic regime 0, 0, 0] , and significant veloc- 
t autocorrelations appear in noncolloidal suspensions 
However, in the quasistatic motion of the present 
experiment, there is no ballistic behavior because grains 
are constantly in contact with each other. Moreover, as 
seen in the velocity autocorrelation shown in Fig. [5] the 
velocities quickly become uncorrelated. The time for the 
correlation to reach zero corresponds to a mean relative 
grain displacement of 0.25 d and occurs within 4 sec- 
onds for this data. Therefore, diffusivities measured for 
5s < t < 30s can be expected to be beyond the corre- 
lated regime and are at times before significant Taylor 
dispersion effects are observed. We show below, using a 
random walk simulation, that the apparent subdiffusiv- 
ity of the radial component in Fig. 0] is due to the radial 
gradient in shear rate, an effect that is not included in 
Eg s . 12141 above. 

Fig. El shows ((RA8) 2 ) and ((Ar) 2 ) for particles ini- 
tially at various distances from the shearing surface. It 
is evident that the diffusivity is larger close to the shear- 
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d<r< 1.5d 
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4.5d < r < 5d 
6d < r < 6.5d 
7d<r<7.5d 
8d<r<8.5d 




d<r< 1.5d 
3d<r<3.5d 
4.5d < r < 5d 
6d<r<6.5d 
7d<r<7.5d 
8d<r<8.5d 




FIG. 6: Mean displacement squared vs. time at different 
distances from shearing wheel. Radial displacements shown 
on top and tangential displacements on bottom. / = lmHz. 



ing surface, i.e. in regions of large shear, as observed 
previously. In addition, the maximum diffusivity occurs 
at r « 2d. The fact that the maximum diffusivity does 
not occur at r — is due to non-diffusive motion of par- 
ticles in contact with the shearing wheel. Since most of 
these particles (r < 2d) are dragged by the wheel at the 
same speed and the mean velocity has been subtracted, 
the fluctuations are smaller. In Fig.0 we show the RMS 
displacements versus time on a log-log scale. The solid 
line shows the expected slope for diffusive behavior. 




Time (s) 




Time (s) 

FIG. 7: Log-log plot of displacement squared versus time for 
data in Fig. |S| Radial displacements given in plot on top and 
tangential displacements on bottom. The straight solid line 
gives the expected behavior for diffusive motion. 



confined to the shear band. 



As seen in Fig. 03 the diffusivity is proportional to the 
local shear rate 7. For this figure, we use the local shear 
rate determined from the slope of Fig. and a diffu- 
sivity that is half the slope of Ar 2 vs. t for t < 30 s. 
The decrease in diffusivity at large shearing rate (i.e. 
close to the shearing surface) is due to particles being 
dragged by the shearing wheel and hence exhibiting bal- 
listic behavior. For the radial and tangential directions, 
D « 0.1 - 0.2 g? 2 7 ks 0.4 - 0.8a 2 7 where a is the particle 
radius, i.e. the scale of the diffusivity is approximately 
given by a 2, y. As a consequence of the exponential tail 
of the velocity profile, the diffusivity also decays roughly 
exponentially, such that the diffusive motion is effectively 



Fig. shows results for diffusivity versus local shear 
rate (r > 2d) for three different rotation frequencies of 
the shearing wheel. The diffusivity is approximately pro- 
portional to local shear rate over a large range of shearing 
rates from separate experimental runs, and Dgg/D rr w 2. 
In Fig. 03, the diffusivity at each data point was di- 
vided by the local shear rate. The resulting data (D/j) 
is roughly constant over 3 orders of magnitude of shear 
rate. The lines show fits for tangential diffusivities Dgg 
= 0.223 7 and radial diffusivities D rr = 0.108 7. The 
data for Dgg is noisier than that for D rr since due to 
the mean flow the magnitude of the tangential motion is 
much larger than radial motion. 
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FIG. 8: Diffusivities D vs. local shear rate j (/ = 1 mHz). 
Diffusivities are proportional to the local shear rate. Close 
to the shearing wheel, at larger shear rates, the diffusivity 
decreases because particles at r < 2 are typically dragged 
continuously by the shearing wheel. The solid line shows D = 
0.200 7 and the dashed line shows D — 0.108 7. 
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FIG. 9: (A) Diffusivities vs. local shear rate for grains at 
r > 2 for three different rotation rates of the shearing wheel. 
In (B), the diffusivities are rescaled by the local shear rate 
7 and d 2 . The lines show fits for Dgg = 0.2237 (solid) and 
D rr — O.IO87 (dashed) which approximately hold for three 
orders of magnitude of shear rate. 



Radial Distance, r (d) 
FIG. 10: Off-diagonal diffusivity D r9 . 

C. Off-diaganol Diffusivity D r g 

The off-diagonal diffusion constant D r g is shown in 
Fig- El This diffusion coefficient is an order of magnitude 
smaller than D rr and Dgg. Away from the shearing sur- 
face, D r g is also negative. This is due to the anisotropic 
force network and will be addressed below. 

IV. RANDOM WALK SIMULATION 

As mentioned earlier, direct comparison of the data to 
Ea. l2l4l is not possible. The fact that v(r) and D(r) decay 
exponentially and the presence of the boundary at r = 
in the experiment are inconsistent with the assumptions 
leading to the moment evolution equations in Section llBl 
An exact solution of the moment equations with the ap- 
propriate boundary conditions and spatial dependence 
of D(r) and v(r) is difficult. However, a numerical sim- 
ulation with appropriate spatial dependence in D and v 
using a random walk model is relatively simple. 

In this section, we present such a simulation in which 
we assume diffusive motion and impose an exponential 
velocity profile and impenetrable inner boundary so as 
to parallel the experiment. To model radial diffusion, a 
walker makes a step each time r with equal probability 
along ±f with a radial step length L r (r) proportional to 
(-^(r)) 1 / 2 , where 7 is the experimentally measured shear 
rate. That is, L r (r) = ci^r)) 1 / 2 , where c\ is a con- 
stant. This imposes a diffusivity D oc L(r) 2 /r cx j(r) 
in agreement with Fig. [5] Radial motion is bounded by 
the shearing wheel, so any step that would move a par- 
ticle through that boundary is automatically forced to 
be a step away from the shearing wheel (i.e. towards 
positive r). Tangential motion is modeled in a similar 
way. At each time step, the walker is advected at the 
experimentally measured mean velocity based on its ra- 
dial position (Fig. 0) and also randomly takes an addi- 
tional step along ±9 with tangential step length Lg(r) 
= C2 A l{r) 1 t 2 This model contains two free parameters, 
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ci and C2, corresponding to the magnitude of D rr and 
Dgg. Equivalently, the fit parameters can be thought of 
as determining an overall scale factor for the data and 
the ratio of radial to tangential diffusivities. Here, we do 
not consider additional anisotropics associated with the 
force chain network. 

Fig. II II shows a mean-square displacement versus time 
for the simulated data. The experimental data from 
Fig. [S] is also included for reference as thin solid lines. 
The simulation is performed assuming Dg/D r = 1.9 {c\ 
= 0.48, C2 = 0.66). As noted, the two free parameters 
set the scale of the radial and tangential displacements. 
We emphasize that the relative magnitudes of the data at 
different distances from the shearing wheel and the ap- 
parent subdiffusive and superdiffusive behavior at longer 
times result from the experimentally measured velocity 
profile. 

Although the initial slope of the lines in Fig.^Jis equal 
to 2_D(r), the long time behavior and deviation from a 
straight line is due to the coupling to the mean flow. 
The horizontal and vertical scales of Figures HJI and ITT1 
are identical in order to compare the long time behav- 
ior. This confirms that the apparent subdiffusion and 
superdiffusion is due to the mean flow. 

In particular, it is clear that the curvature of the D rr 
data for particles close to the shearing surface, which 
appeared to be subdiffusive, arises from the gradient in 
local shear rate. That is, the grains next to the wall dif- 
fuse away to a region of slower shearing rate and, once 
away from the wall, diffuse more slowly. If this gradient 
is removed from the simulation (i.e. 7 is assumed con- 
stant), the lines become straight with approximately the 
same slope. The presence of a wall (r = 0) also tends 
to decrease the diffusivity at small r, but this is a much 
less pronounced effect than that of the gradient in shear 
rate. Note in Fig. ^2 that simulated grains close to the 
shearing wheel show an apparent superdiffusive behavior 
due to Taylor dispersion. 

We note that in the experiment, grains at r < 2 are 
generally dragged by the shearing wheel which tends to 
decrease the apparent radial diffusivity and add to the 
effects of Taylor dispersion (which accounts for the slight 
difference in the magnitudes of Fig. l6land lllfl . However, 
Fig. ^3 reveals that the dominant effects are the shear 
gradient and Taylor dispersion. Thus, the main features 
of the apparent subdiffusive and superdiffusive behavior 
are observed even though the simulation does not include 
the effect of ballistic motion due to dragging of particles 
by the shearing wheel. 

We show the measured diffusivities for the experiment 
and the random walk simulation together in Fig. 1121 
There is very good agreement except for r < 2, where 
the simulation overestimates the diffusivity. 




Time (s) 

FIG. 11: Mean displacement squared vs. time for simula- 
tion of random walk in which the experimentally measured 
velocity profile is imposed. Radial displacements shown on 
top and tangential on bottom. The thin solid lines show the 
experimental data from Fig. |5] 
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FIG. 12: Diffusivity vs. radial distance from the shearing 
surface for experimental data and the random walker model. 
The deviation between simulation and experiment for r < 2 is 
due to particles being dragged nondiffusively by the shearing 
wheel. 





FIG. 13: (A) Force chains imaged using photoelastic grains. 
The shearing wheel is marked by the white line and is rotating 
to the left. (B) Image (A) is rescaled from polar to cartesian 
coordinates such that the shearing wheel is at y = 0. (C) A 
2D autocorrelation of (B) which characterizes the orientation 
of the strong force network. 



V. EFFECTS DUE TO ANISOTROPIC FORCE 
NETWORK 

In the previous discussion, we tacitly assumed that the 
natural coordinate sytem for diffusion measurements is 
set by the radial and tangential directions, correspond- 
ing to the anisotropy of the imposed shear. However, 
dense systems, unlike dilute rapid flows, have anisotropic 
force networks due to imposed shear which are in general 
at a different orientation from the flow direction. This 
is seen by using the photoelasticity of the grains to im- 
age the force chains, as in Fig. 1131^. This figure shows 
a typical case where the force chains are preferentially 
oriented to oppose the motion of the shearing wheel, at 
an angle that is intermediate between the f and 9 di- 
rections. We might expect the minimum diffusivity to 
be oriented along the mean direction of the force chains. 
To determine this angle, images such as Fig. can 
be transformed from polar to cartesian coordinates, such 
that the shearing wheel is located at y = (Fig.^j3). In 
this figure, each vertical line corresponds to a radial line 
in the original image. Since the curvature of the wheel is 
relatively small, the transformation is not dramatic and 
distortion of the image is small. A 2D autocorrelation of 
image IT3T 3. then provides a measure of the mean force 
chain orientation 1)130). The mean angle <f> of the force 
chains fluctuates strongly in time around a mean value 
of 20 - 30° relative to f. 

To determine the angular dependence of the diffusive 
motion, we locally project displacements at each time 
step onto an axis rotated at an angle <fi from the radial 
direction, as sketched in Fig. ^] That is, for each step, 
the displacement is locally parameterized in terms of ra- 
dial and tangential components (relative to the center of 
the shearing wheel) and then are locally projected onto 
an axis rotated by <f> relative to the r direction. We then 
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FIG. 14: We project displacements at each time onto an 
axis rotated at an arbitrary angle </> degrees from the radial 
direction. is measured versus <j> and plotted with the 

minimum diffusivity indicated by the filled circles. 
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FIG. 15: The angle cj> corresponding to the minimum dif- 
fusivity 4> m in is measured from Fig. 1141 The increase with 
radial distance corresponds to the minimum diffusivity be- 
coming more aligned with the mean direction of force chains. 



use the <p components of the trajectories to measure a 
diffusivity, D^, along this direction. On the right side 
of Fig. El we show versus 4> for grains at different 
radial distance r. Again, the diffusivities decrease with 
distance from the shearing wheel. In addition, the direc- 
tion of minimum diffusivity 4> m in, marked by the solid 
circles, changes with distance from the shearing wheel. 

In Fig. 1151 we show the angle of minimum diffusiv- 
ity <pmin versus distance from the shearing wheel which 
we determine by fitting a parabola to in the re- 

gion <t) m i n ± 30°. Close to the shearing surface, at high 
shear rates, the minimum diffusivity is in the radial di- 
rection, corresponding to the minimum expected based 
on the imposed shear direction. At larger distances r, 
the minimum shifts towards the direction of the mean 
force chain orientation. In other words, outside of the 
immediate vicinity of the shearing wheel, the anisotropic 
force network affects particle motion and must be taken 
into account in order to properly describe the diffusive 
motion in dense granular systems. 

We use the same procedure for deducing the angular 
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FIG. 16: The angle of minimum diffusivity cf> min VS. V IH611- 
sured according to Fig. 1141 for multiple data sets. Again, the 
orientation of minimum diffusivity shifts towards a direction 
that corresponds to the mean force chain direction. The same 
analysis is performed on the random walk data (o and linear 
fit) which does not model the anisotropic force network. 
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FIG. 17: Diffusivities determined from velocity autocorre- 
lations using Eq. The shearing wheel frequency is / = 1 
mHz. The lines show fits to the tangential (solid) and radial 
(dashed) diffusivities in Fig. I12l for comparison. 



dependence of the diffusivity on a number of independent 
data sets, including the data used to create the velocity 
profile of (Fig. |3J) and on two additional data sets with 
different shearing rates. We show data in Fig. 1161 only 
for velocities that are properly resolved, as in Fig. [3] 
Although there is some significant variability from one 
data set to the next, there is a clear trend in which <j) m i n 
shifts towards the direction of force chain orientation as 
r increases. When we perform the same analysis on the 
simulated data, in which there is no force network, (j) m in 
does not increase above 0° as indicated by the open cir- 
cles. The fact that these points are negative is addressed 
in the following section. 
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VI. DIFFUSIVITY AND VELOCITY 
AUTOCORRELATION FUNCTIONS 

Diffusivities {D rr and Dgg) can also be determined us- 
ing velocity autocorrelations from the expression 



D xx 



{v x {t)v x (t + r))dr, 



(7) 



where the velocity at each time step is simply defined 
as v x (t) = (x(t) — x(t — 1))/ At. The integral must be 
taken over times long enough to extend beyond the initial 
correlated region. Thus, for this data, 



D , 



At 



N 

£< 

dt=0 



v x (t)v x (t + dt)}, 



(8) 



where we use a cumulative sum of the autocorrelation 
(Fig. [SJ. After an initial transient of about 5 s, the curve 
fluctuates around a constant value. The average of the 
data well after the transient (85 - 185 s, for this data set) 



FIG. 18: Off-diagonal diffusivity determined by integrating 
the velocity cross-correlation. Data determined from the slope 
of displacement squared versus time plots (Fig. 1101 is shown 
for comparison. D r g for the random walk simulation is also 
shown. 



is taken as the diffusivity. Fig. El shows the diffusivities 
determined from the velocity correlations. For the most 
part, they agree quite well with diffusivities determined 
from the displacement squared versus time data which 
are indicated by lines showing fits to the data of Fig. IT21 
The exception is for tangential diffusivities of particles 
adjacent to the shearing wheel. This discrepency is due 
to the fact that that there is some non-negligible correla- 
tion in the tangential velocity after t — 30s for particles 
in contact with the shearing wheel. 

The off-diagonal diffusivity must be determined using 



D 



V = 2 / + (*')»<**' (9) 
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which for large t can be rewritten as L8jj 



Di 



{vi (t)vj(t + r))dr, 



(10) 



assuming the motion is statistically stationary over time. 
This reduces to Eq. [7| for the diagonal terms, i = j. The 
off-diagonal diffusivity, D r g is shown in Fig. llSl along with 
data determined using the displacement squared versus 
time plots (Tig. ITU)) . Since the magnitude of D r g is small, 
the fluctuations lead to larger noise in these results. 

Further from the wheel, where the minimum diffusivity 
shifts to larger angles, the cross-correlation term is nega- 
tive. This indicates that motion along is anticorrelated 
with motion along f , which agrees with a decrease in dif- 
fusivity along positive (f>, and a shift in minimum diffusion 
angle towards positive tj>. To emphasize the effect of the 
anisotropic force network, we contrast the experimental 
results with data from the simulation, where there is no 
force network effect and D r g is always positive (Fig. I18|l . 

The fact that D xy is positive in the simulation is due 
to the velocity gradient, an effect that was also observed 
in previous measurements of the cross-term 0, 0, . 
Note that a positive D xy in the present data corresponds 
to a negative D xy in previous results. This difference in 
sign is due to the fact that the authors of [l7lll^ . l26l | used 
the convention that v x increases with y. By contrast, we 
have chosen the natural experimental coordinate sytem 
y = r such that y = at the boundary (shearing wheel) 
and increases into the bulk, i.e. v x decreases with y. 

Returning to the present studies, the fact that D xy is 
positive for the simulation is the origin of the slightly 
negative (j)min seen in Fig. 1161 We emphasize that the 
positive (fimin and negative D xy in the experimental data 
have the opposite sign from the simulation. This differ- 
ence in sign is due to the presence of the strong force 
network in the experiments, an effect that is absent in 
the simulations. 



VII. LEVY FLIGHTS AND VORTICES 

We also observe examples of correlated motion and 
trajectories similar to Levy nights [35| which could con- 
tribute to non-diffusive motion. In fact, the dense 2D 
packing leads to caging and coordinated motion such that 
neighboring grains tend to move together. This differs 
from more dilated flows in which collisions are the source 
of fluctuating motion. The fact that the present system 
behaves diffusively on average indicates that long-range 
correlated motion is sufficiently rare and random over 
time that the mean behavior is not affected. 

It is interesting to ask whether such novel behavior as 
Levy flights occur in our system and whether they are 
important. Levy flights are random walks in which oc- 
casional large steps, or flights, are observed, such that 
apparent Brownian motion on smaller scales is punctu- 
ated by large displacements. They also have the property 
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FIG. 19: (A) Trajectories in the shear band are shown for 
a fixed time (1000 s). (B) A few individual trajectories from 
the dashed box in (A) are shown. Motion similar to Levy 
flights is occasionally observed (r ~ 4 — 5d), in the region in- 
dicated by the arrow. Small fluctuations of particle position 
are observed with occasional larger scale advection. (C) Tan- 
gential displacements versus time for 7 individual trajectories 
(r « 2.5, 3, 3.5, 4, 4.5, 5.5, 6.5d) chosen to highlight L'evy-like 
motion. Trajectories closer to the shearing wheel are dis- 
placed upwards at t — for clarity. 



that the variance of the step size (L 2 ) and therefore the 
diffusivity (oc L 2 /t) are infinite. This situation can be 
realized if the probability of the walker making a step L 
is given by a power-law P(L) oc L~ a where 2 < a < 3. 
This is in contrast with gaussian or exponential distribu- 
tions of L in which case large steps are much more rare 
and the variance is finite. 

In Fig. we show trajectories for particles in the 

shear band. Each line shows the trajectory of a single 
particle over the same time period (1000 s). Next to the 
shearing surface (bottom of image), particles travel rel- 
atively fast compared to particles outside of the shear 
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band (top of image), which fluctuate around effectively 
stationary positions. Fig. 119b shows a few particular tra- 
jectories from the dashed region of Fig.lT^K. At the edge 
of this band (r « 4 — 5d), we see trajectories that are 
reminiscent of Levy flights, in which relatively large dis- 
placements occur between periods of fluctuating motion 
on a smaller scale. To observe this better, in Fig.lliJC. we 
plot tangential displacement versus time for seven trajec- 
tories at different distances from the shearing wheel. The 
data has been smoothed over a 5 s running window. We 
note that these are not necessarily typical trajectories, 
but have been chosen to elucidate the presence of Levy- 
like behavior. In this plot, regions of small fluctuations 
(nearly flat lines) are separated by faster jumps along the 
mean flow direction. 

These data can be compared with an example from a 
rotating flow fluid experiment |3rl e.g. Fig. 7] from which 
similar data were obtained. Although we observe some- 
what similar motion, the trajectories in Fig. 1191 have a 
different origin. In the fluid case, particles exhibit flights 
between periods in which they are trapped by vortices. 
In the granular system, grains become trapped as they 
move farther from the shearing surface and remain effec- 
tively trapped until they move closer to the wheel. In 
addition, it is common for local rearrangements involv- 
ing 10-20 grains to occur intermittently. (A similar effect 
may also account for the Levy distributions of trapping 
times for observations of grains deposited on sand piles 
|37|-1 With the present data, e.g. Fig. 1191 we do not have 
sufficient statistics to determine whether the trajectories 
exhibit Levy scaling because flight-like trajectories are 
rare. 

We also occasionally observe cooperative motion as in 
Fig. I^nj which displays particle trajectories over a 25 s 
window in which the greyscale level indicates time (light 
grey = early time, dark = later time). In the lower right, 
there is a region of locally correlated motion. A tran- 
sient vortex is present in the upper left. Although corre- 
lated motion is common, since motion in a dense packing 
requires motion of neighboring grains, vortices are rare 
events. In addition, unlike vortices in fluids, there are no 
inertial effects and granular vortices appear to quickly 
dissipate without affecting the long time behavior of the 
grains. 

VIII. CONCLUSIONS 

To conclude, we find that granular motion in dense 
shear flows is diffusive with a self-diffusivity proportional 
to the local shear rate (D ~ 7a 2 where a is the particle ra- 
dius). However, the diffusion tensor D is anisotropic due 
to underlying anisotropies in both the velocity field and 
force network. The velocity anisotropy leads to a tangen- 
tial diffusivity that is about double the radial diffusivity. 



The anisotropic force network dominates the local diffu- 
sivity outside the immediate vicinity of the shearing sur- 
face, and leads to a minimum diffusivity approximately 




FIG. 20: Granular vortices are occasionally observed in plots 
of particle trajectories over 25 s. Particle position is indicated 
by a dot with darker grayscale levels used for later times. 

along the mean force chain direction. This latter fea- 
ture has not been observed in more rapid flows, to our 
knowledge, and is a property of dense granular systems. 

Motion can appear subdiffusive due to the decreasing 
shear rate away from the shearing surface or superdiffu- 
sive due to Taylor dispersion effects. A simple random 
walk model which reproduces the apparent anomalous 
diffusion indicates that the underlying motion is diffu- 
sive. Using the experimentally measured velocity profile, 
assuming Dg/D r w 2, and choosing an overall multi- 
plicative scale factor, the simulation closely matches the 
experiment, including the long time behavior which is 
affected by the gradient in shear rate and Taylor disper- 
sion. The simulation also highlights the effects of the 
anisotropic force network. Differences in the sign of D xy 
between simulation and experiment are associated with 
the anisotropic force network. The same is true for the 
orientation of the minimum diffusivity D^. Velocity au- 
tocorrelation plots show that motion in dense granular 
flows quickly becomes uncorrelated and there is not a 
distinguishable ballistic regime before diffusive behavior 
dominates. 

Examples of correlated motion, such as vortices, and 
trajectories similar to Levy flights are also observed. 
However, these effects are sufficiently intermittent and 
random that the system behaves diffusively. 
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